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ABSTRACT 

Recent detections of GeV photons in a few gamma-ray bursts (hereafter GRBs) by 
Fermi-hAT imply huge bulk Lorentz factors to avoid a large 77 optical depth at high 
energy. Estimates can be as high as F ~ 10^ in the most extreme cases. This puts severe 
constraints on models of the central engine and the jet acceleration in GRBs. These 
estimates are however obtained from a simplified single zone model. We present here 
a more realistic calculation of the 77 opacity which takes into account the time, space 
and direction dependent photon field existing in an outfiow with several relativistically 
moving emitting zones. The formalism is very general and can be applied to many 
models of the prompt GRB emission. We present results obtained for a numerical 
implementation in the framework of the internal shock model. We show that (i) the 
minimum Lorentz factor Fmin in bright Fermi-LAT GRBs is reduced by a factor 
~ 2 — 3 compared to previous estimates if the GeV and MeV emission are produced in 
the same region, and by an additional factor ~ 2 — 8 if the GeV emission is produced 
at larger radii. We provide an improved approximate formula for Fmin which is in good 
agreement with our numerical results and which can be directly applied to Fermi-LAT 
GRB data; (ii) a delayed onset of the GeV emission can be due to the time evolution of 
the opacity in a GRB outfiow. As an illustration of these first two results, we present 
a synthetic GRB that reproduces most features of GRB 080916C with a mean Lorentz 
factor of ~ 340, an optically thin regime for 77 opacity at 3 GeV in time bin 'b', a 
variability timescale of ~ 0.5 s in the MeV lightcurve and a delayed onset of ~ 5 s of 
the GeV emission; (iii) the 77 opacity can smooth the short timescale variability in the 
GeV lightcurve. This last result implies that the observed variability at high energy is 
not necessarily a good test to distinguish between an internal and an external origin 
for the GeV emission in GRBs. 

Key words: gamma-ray burst: general; gamma-ray burst: individual: GRB 080916C; 
radiative transfer; radiation mechanisms; non-thermal. 



1 INTRODUCTION 

The combination of a short timescale variability with a huge 
radiated gamma-ray energies in gamma-ray bursts leads to 
the well-known compactness problem: 7-ray photons should 
not escape due to 77 annihilation. This contradiction is 
solved by assuming that the e mitting region has an ultra- 
relativistic motion (|Reeslll966l ). This leads to an important 
constraint on the minimum Lorentz factor of the outflow in 
cases where the GRB spectrum do es not show any evidence 
for a high-energy cut ofi^ (see e.g. IBaring fc H arding 19971: 
iLithwick fc Sarill200ll '). or to a direct measurement of the 
Lorentz factor if such a cutofi^ is iden tified, as possibly in 
GRB 090926 A l|Ackermann et al.ll201ll ). 
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Since the launch of Fermi in June 2008, the LAT 
instrument has detected high energy photons above 10 
GeV in a few GRBs. The observed 7-ray spectrum often 
remains consistent with a Band function covering the GBM 
and LAT spectral ranges without any evidence of a high 
energy cutoff which could be identified as a signature of 77 
annihilation (with poss ibly the exception of GRB 090926A, 
lAckermann et al.l I2OIII ). This extension by Fermi of the 
observed spectral range upper bound from a few MeV to 
a few tens of GeV implies constraints on Fmin which are 
mu ch more severe than the ones obtained previously (see 
e.g. iRacusin et al.ll201ll V In a few cases Fmin has been esti- 
mated to be of the order of 1000. For instance, in the four 
brightest Fermi-LAT GRBs, the estimated minimum values 
are F^in ■r^ 887 in G RB 080916C jAbdo et al.ll2009bl ). ~ 10^ 
in GRB 090902B (|Abdo et al.l l2009al ). ~ 1200 in GRB 
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090510 llAckermannet aLll2010l )). ~ 720 in GRB 090926A 



IjAckermann et alj 12011 



These extreme values put severe 
constraints on the physics of the central engine which should 
be able to strongly limit the baryon load in the outflow. 
However these minimum Lorentz factors were obtained 
from a simplified "single zone" model where the spatial and 
temporal dependencies of the radiation field are averaged 
out. This calls for more realistic calculations of the 77 opac- 
ity in GRB outflows. We present here a detailed approach 
taking into account the exact time, space and direction 
dependent radiation field in outfiows with several rela- 
tivistically moving emitting regions, such as internal shocks. 

The paper is organized as follows. In Sect. [2l we de- 
rive the exact formula for the 77 opacity in complex and 
variable GRB outflows, we describe its dependency on sev- 
eral physical parameters such as the emission radius or the 
Lorentz factor and we compare the exact calculation to the 
approximate formula which is broadly used in the literature. 
Finally we validate our approach by a comp arison with the 
semi-analytical study of lGranot et al.l l|2008l ). Sect.[3]is de- 
voted to the investigation of several consequences of 77 an- 
nihilation in GRBs in the framework of the internal shock 
model. We discuss successively the constraints on the min- 
imum Lorentz factor, with a focus on GRB 080916C, the 
spectral shape of the high-energy cutoff, the delayed on- 
set of the GeV emission and the smoothing of the short 
timescale variability in GeV lightcurves. We compare also 
in this section the 77 opacity with other expected sources of 
opacities. Sect. |4] presents a discussion of models where the 
GeV emission is produced at a larger radius than the MeV 
main component. Finally, Sect. [5] is the conclusion. 
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Figure 1. Geometry of the 77 interaction of a high energy 
photon with a spherical flash. All quantities are defined in 
the source frame. A high energy photon is emitted at radius _Re 
and time te (point E) in a direction ©e relative to the radial 
direction, and interacts with a flash emitted at radius Rq and 
time to- The figure describes the interaction geometry at a given 
time tj > to : F is the emission point of the low-energy seed photon 
and I the interaction point. The angles 8, (f> and 5 are respectively 
the emission colatitude, longitude and Doppler angle of the seed 
photon. Note that the path EI (high energy photon) is contained 
in the plane of the figure, which is not the case for the path FI 
(seed photon). The radius of the interaction point I is R[ and, 
from I, the arrival direction of the seed photon and of the HE 
energy photon make respectively an angle oj and O with the 
radial direction. The distances from F to I and from E to I are 
respectively s and £, leading to ti = te +i/c = io + s/c. Finally ^ 
is the angle between the directions of the two interacting photons 
(hereafter the interaction angle). Note that, at a given time tj, 
all the seed photons come from the same emission colatitude 0, 
leading to the same values of oi, s and S = 6 + ai. However, the 
angle ip also depends on </>. 



2 CALCULATION OF THE 77 OPTICAL 
DEPTH 

2.1 General 77 opacity formula. 

The 77 opacity r-^^ seen by a photon of high energy -Ehe 
propagating in a relativistic outflow is given by: 

r^-y(£'HE) =[d£ IdQ. [ dE nn{E)a^-,{EHE, E,iP) (1 - 



and y is defined by 



■sip) 

(1) 
All the physical quantities are measured in the fixed frame 
associated to the central source ("laboratory frame", here- 
after the source frame): E is the energy of the interacting 
field photon, -0 represents the interaction angle between the 
HE photon and the interacting photon. The energy thresh- 
old Ec of the field photon above which 77 annihilation can 
occur is 
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-Ehe (1 — cost/;) 
and the 77 cross section a-fj is given by 

a^-yiy) = aTg{y), 
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for E ^ Ec a.nd y = ior E ^ Ec. Finally nn{E) is the 
photon field distribution [ph • cm~^ ■ keV 
location and time. 



""^ ■ sr ^] at a given 



Eq.((T} is made of a triple integral: the d^-integration is 
done over the path of the GeV photon from its emission loca- 
tion to the observer, the dn-integration is done over the solid 
angle distribution of the interacting photon field surround- 
ing the GeV photon whereas the d_B-integration is done over 
its energy distribution. Eq.([T| is general and can be applied 
to any photon emitted at a given location and time with a 
given direction of propagation in the GRB outflow. 

2.2 77 opacity created by a "spherical flash". 

2.2.1 Calculation of the 'yy opacity 

The kernel of our study is the calculation of the 77 opacity 
created by a spherical flash. Using this kernel, it is possible 
to model the case of a propagating radiating spherical front 
(representing for example a shock wave) by the succession 
of many spherical flashes. One critical step is the exact cal- 
culation of the photon density nn (see Appendix \K\ . The 
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emissivity [erg • cm~^ ■ s~^ • Hz~^ ■ sr~^] in the comoving 
framqj of a spherical flash occurring at a radius Ro and 
time to (source frame) is given by : 

^'•^' = JZJ^^-^:!^^(4-)5{r~R.)5{t^to), (6) 
47r47ri?(, To u^f^ \'^p,o J 

where f^ad is the total radiated energy, Fq — l/-\/l — /Jq is 
the Lorentz factor of the material moving with velocity /3oc 
and emitting the spherical flash, v' is the photon frequency 
in the comoving frame, z^p g is the comoving peak frequency 
of the spectrum, and B is the normalized spectral shape 

a;°B{x)dx = i). 

The geometry associated to a spherical flash is described 
in Fig. [T] and the calculation of the created photon field 
nn{E) is given in Appendix [A] The resulting general ex- 
pression of the 77 opacity created by a spherical flash is 
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where _Ep q = hu'^Q is the comoving peak energy, i/; is the 
interaction angle given by 

cos ^j) = cos O cos ai + sin O sin ai cos (j) , (8) 

D is the Doppler factor 

1 



V = 



Fo(l-/3ocos5) ' 
and the dimensionless constant tq is defined by 

CT^rad 
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If the high energy photon is emitted at radius Re and 
time tc in a direction 0c relative to the radial direction (see 
Fig- Hi the quantities appearing in Eq.Q are given as a 
function of the length I traveled by the high energy photon 
from its emission point by 



RI = rI + l^ + 2RdcosQe, 
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Physical quantities with a prime are measured in the comoving 
frame. 



In Eq.lO, the d^- integration is limited to the range of 
allowing an interaction, i.e. fulfilling the condition 



\Ro-Ri\ < s^i?o + i?i. 



(17) 



— Case of a photon emitted on axis. If the high energy 
photon is emitted on axis (Oo = = and ■(/; = a/) the 
axial symmetry gives a simplified formula : 
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— Case of a flash having a power-law spectrum. If 

the spectrum of the flash emission can be approximated by 
a single power-law of photon index /3, the spectral shape 
function has the following form B{x) — Aqx^'^^ for x ^ 1 
with Ao = -(2 + P), and B{x) = for a; < 1. Then Eq.Q 
can be further simplified: 
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where I(/3) is a finite number for /? < : 

^(P)-l \iJy2)2+, 9iy)dy- (22) 

For /3 — —2, —2.2, —2.3, —2.5 and —3, this integral equals 
respectively I(/3) ~ 0.091, 0.078, 0.072, 0.063 and 0.046. 

— Case of a photon emitted on axis interacting with 
a flash having a power-law spectrum. Finally if the two 
simplifications are combined, the r-y^ formula becomes: 
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(23) 



2.2.2 Approximate formula 

Most estimates of the minimum Lorentz factor in GRB 
outflows are based on approximate formulae corresponding 
to the case of a high-energy photon emitted on axis and 
interacting with an isotropic radiation field over a length 
£' ~ -R/Fo in the comoving frame of the emitting material 
assu med to move with a constant L orentz factor Fq (see e.g. 
i Lithwick fc Saril l200l[: iGupta fc Z hang 200^; lAbdo et al.l 
l2009bl : IZhang fc Pe'eij|2009l ). In manv GRB models, such as 
internal shocks, a more realistic approximation corresponds 
to the situation where the seed photons are produced by 
a flash emitted at radius Ro by a single thin shell moving 
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with constant velocity Po, and where the high-energy pho- 
ton is emitted on-axis at -Re and interacts with the resulting 
anisotropic radiation field. Then, the law of motion of the 
thin shell is Re — Ro ~ l3oc{ta — to) and the high energy 
photon can interact with the seed photons from the fiash 
only if 



-Ro < ^o < 



2-Ro 



4rgi?o . 



(24) 



In this case, an approximate formula for the 77 opacity given 
by Ea. (|23p can be obtained (see Appendix [B| : 
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(25) 

To compare with observations, the radius _Ro is usually es- 
timated by FocAfvar, where Aivar is the observed variabil- 
ity timescale, and the isotropic equivalent radiated energy 
frad is deduced from the observed photon fluence T{E) = 
To {E/Ep^oY [ph • cm-2 • keV^], where E^^o = ^oE'^^ is 
the observed peak energy. This leads to 



used bv lAbdoet"ai] (|2009bh differs from ATf*^ by less than 
0.5%. The coefficient 

used by iLithwick fc Saril (|200ll ) is less accurate, with a 
difference by a factor 1.5 to 2 in the range — 3 < /? < — 2. 

For a typical value /3 = —2.3, the differences between 
the two assumptions (isotropic radiation field vs anisotropic 
field created by a flash) is noticeable. We have A'^°®/A"o — 
14.7 and K^^°^ / Ko ~ 7.3, already leading to some differ- 
ences in the estimate of the minimum Lorentz factor in 
GRB outflows. Note that in all these single zone models, 
the opacity is built over a length I ^ Rm the source frame 
(see Appendix |B] and see also ^3.3. l|l , so that the origin of 
the difference is mainly the geometry of the radiation field 
(isotropy vs anisotropy in the comoving frame). However, as 
described in i]3.3l the estimate of Fmin is even more affected 
by dynamical effects, that are not included in these single 
zone models. 
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where D is the distance of the sourcqj Finally the 77 opacity 
reads 

X 2 / / 2\2 \ 1+/3 

Tj^f ~ K (Re) ar \ —r- — Ep,oToTq 
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where 
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(28) 



is a coefficient that depends on the location where the 
high-energy photon is emitted. Most previous studies con- 
sidered usually the limit Re — ^ Ro (all photons are radi- 
ated in the same region). This leads here to a coefficient 
Ko ^ K (Re -> Ro) given by 



A"o 



(29) 



Compared to the simple model assuming an isotropic radia- 
tion field in the comoving frame of the outfiow, the approxi- 
mate formula derived here for an anisotropic field produced 
by an impulsive flash shows the same scaling with the pa- 
rameters of the problem, but a different normalization Ko- 
ISvenssonI (jlOSTJ) has derived the exact coefficient for the sin- 
gle zone isotropic model and gives the following accurate 
approximation (error less than 0.3% for — 7 < /? < — 1) 



2. 2., 3 A case study 

To better understand the different effects at stake when 
studying the 77 opacity in relativistic outflows, it is worth to 
focus first on the case of a single spherical flash for which we 
have obtained an approximate estimate in H2.2.2I Let us con- 
sider a spherical shell expanding at a constant Lorentz factor 
Fo that emits a photon flash at a radius -Ro = (Fo/100)'^ 10^* 
cm. The total radiated energy (source frame) is £rad ~ 2-10^^ 
erg and the spectrum is a broken power-law with a comoving 
peak energy E'p^ = (1 MeV) /Fo, a low energy photon index 
a = —1, and a high energy photon index /3 = —2.5. With 
these assumptions, the observed pulse has a fixed duration 
Atobs ~ Ro/{2Toc) ~ 0.17 s, a fixed bolometric isotropic 
equivalent energy f^ad and a fixed observed peak energy 
~ Fo-Bp — 1 MeV. In addition, we assume that the same 
shell, expanding at a constant Lorentz factor Fq, is emitting 
a high energy photon at a radius Re ^ Ro and a colati- 
tude 0c. Eq.([7]) allows to compute the optical depth for 77 
annihilation t-,^ seen by this high energy photon, as a func- 
tion of its energy -Ehe- This photon is detected by a distant 
observer at time t™ given by 

,e .., .flash I ^c JlO /., r^ \ ^c /oo\ 

iobs - iobs H ^7^2 ^^ (1 ~ COS Bo) — , (33) 

where tob? is the time at which the first photon from the 
fiash is observed. 

The four panels of Fig. [2] give an overview of the obtained 
results showing the different dependencies of r-yj. 
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2 We do not include here the corrections due to cosmological 
redshift. 



Dependence on the emission radius. The top left panel 
of Fig. [2] gives the evolution of r^-y as a function of Re for a 
photon emitted on axis (Oe = 0). The opacity r^-y remains 
constant for 






''obs 



iflash 
^obs 



Afo 



<0.2 



(34) 



and then decreases as Tjj oc Re , in excellent agree- 

ment with the approximate formula given by Ea. (|27p . The 
top right panel of Fig. [2] shows the same curves T-,^{Re), 
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Figure 2. Opacity created by a spherical flash. In each panel, the 77 opacity is plotted as a function of the high energy photon 
emission radius. In the first three panels, the Lorentz factor is To = 100. Upper left panel: opacity seen by a 1 GeV photon emitted on 
axis. The numerical calculation is plotted in solid line and the approximate formula given by Ea. ll25l l in dashed line. Upper right panel- 
opacity seen by a 1 GeV photon for different emission colatitudes SeTo = (black, solid line); 0.03 (black, dashed line); 0.1 (black, 
dashed dotted line); 0.3 (black, dotted line); 1 (gray, solid line); 3 (gray, dashed line); 10 (gray, dashed dotted line). Bottom left: opacity 
seen by a high energy photon emitted on axis for EuE = 1; 10; 100 MeV; 1; 10; 100 GeV ; 1; 10; 100 ; 1000 TeV. Bottom right: opacity 
seen by 1 GeV photons emitted on axis for different Lorentz factors Fq = 100; 200; 300; 400. 



but for different emission colatitudes ©e. The evolution is 
more complicated but remains dominated by a transition 
from a plateau when Ro is still close to Ro, to a steep decay 
at large radii. As explained in Appendix [B] the plateau is 
due to a favorable configuration when R^ — >■ Rq: both the 
Doppler angle 5 and the interaction angle 1/) are close to 
l/Fo at the beginning of the propagation of the high-energy 
photon; at larger emission radii Re, the Doppler angle 
increases, leading to a reduction of the specific intensity, 
and the interaction angle decreases, leading to a higher 



energy threshold for the seed photons and a lower effective 
cross section (1 — cos-i/)) a^^. 

Dependence on the emission latitude. As seen in the 
top right panel of Fig. [21 the value of the opacity in the 
plateau discussed in the previous paragraph depends on the 
emission colatitude Q^ of the high energy photon. When Q^ 
is increasing, the dominant effect is initially an increase of 
the interaction angle ^ at the beginning of the propagation, 
which leads to an increase of the opacity. However, when 
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the colatitude Qe increases even more, this effect is reduced 
by an increasing Doppler angle : the opacity reaches a 
maximum for Qa — 0.3/To and then decreases. At very high 
colatitudes, the opacity is strongly reduced compared to 
the on-axis photon: e.g. by a factor ~ 30 for Oc — IO/Fq. 

Dependence on the photon energy. The bottom left 
panel of Fig. [2] gives the evolution of r-y^ as a function of Re 
for different photon energies Ehe- The curves are similar to 
the on-axis case discussed above. However the value of r-y^, 
in the initial plateau evolve with -Ehe. For ^ ~ 1/Fo (low 
emission radius Re), the threshold for 77 annihilation is 
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(35) 
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Then, for Ehe < (2Fomoc'') V-Ep,o =i lOGeV, the thresh- 
old energy Ec remains in the high energy branch of the 
spectrum and the dependency is simply r-,^ ex E^^ 
with —(1 + j3) = 1.5 for the present case. Then, for _B > 10 
GeV the threshold E^ begins to enter the low-energy 
branch of the spectrum and the dependency progressively 



evolves from t^i oc E, 



-a+m 



to r-7 



c< EttL " because 



photons contributing mostly to r^^ have an energy close 
to Ec. Due to the particular choice of a — —1, we have 
— (1 + a) = and the 77 opacity r-y-y becomes independent 
of She at very high energies (see the saturation effect ob- 
served in the middle panel of Fig. [2] when J5he is increasing). 

Dependence on the Lorentz factor. The bottom right 
panel of Fig. [2] gives the evolution of r-f-y as a function of 
Re for different Lorentz factors Fq. In this case the scal- 
ing is simply Tjj oc Fq . As discussed in Appendix iBl 
this scaling comes from the combination of three effects: 
(i) the threshold energy goes like Ec oc Fq so that Tj^ 
is modified by a factor Fq' ' (see Ea. (|23p ) (ii) Atobs is 
fixed and Ro oc Fq, which modifies t^j by another fac- 
tor Fq : the specific intensity is diluted by a factor R^^ 
(see Eg. HlOp ) but the opacity is built on a typical length 
I ~ Ro (see Appendix |B| ; (iii) the geometrical factor 
(1 — cosip) ~ I/2F0 brings another factor F~^. Finally, the 



dependency is t-,^ oc Fq' 



2(l+,9)-2-2 



= F, 



2(13-1) 




2.3 Towards a complex photon field geometry 

Using the present approach, it is possible to treat any com- 
plex photon field configuration as long as it can be modeled 
as a collection of spherical photon flashes. The contribution 
from a flash emitted at Ro, to by a expanding sphere with 
Lorentz factor Fq, with a spectrum defined by the comoving 
peak energy _Ep,Q and the shape B and with a total radiated 
energy £^rad, to the 77 opacity seen by a high-energy pho- 
ton of energy -Ehe emitted at Re, te in direction Oe can be 
formally written as 

^r D + r-l \ CTT^rad 

r^7 (iiHE,-Kc,tc,fc)o) = . p2r g/ 



/ d£ J^ [£; Eb_e, Re,te, Qc', Ro, to, Epo,B\ 



where the formal function J^ can be deduced from Eq.Q, 
and Eq.® to Ea. (|17|) . In a more general situation, the 77 



opacity seen by a high energy photon is given by 

r-,, {EnE,Re,tc,e. ^ -^ f-^^ <^TA.ad,.(i) 



E/..; 



4^i??(i)F.(i)i5;,,(f) 

X f d£J^[£; EHE,Re,te,ee; R^{t),t,E'p^,{t),B^{t)] , 

(37) 

where the sum on i is on all emitting regions in the out- 
flow (such as the shocked region behind a propagating shock 
wave), the di-integral is computed over the whole propaga- 
tion of the considered emitting region (t being measured in 
the source frame), and the elementary contribution to the 
opacity of the propagating emitting region i at time t is 
given by the same formal expression as for a spherical flash, 
replacing the radiated energy frad by the radiated luminos- 
ity /Irad,i(i) at time t, and considering all other properties 
(radius, Lorentz factor, emitted spectrum) at time t. Finally 
synthetic light-curves and spectra seen by a distant observer 
can be computed by integrating the flux from each emitting 
region on equal-arrival time surfaces, including at high en- 
ergy the effect of r-y-y . 



2.4 Validation of the model 

Before dealing with more complex dynamical configura- 
tions within the internal shock framework, the validity of 
our model has been tested on a simple single-shock case 
with a comparison to the previous semi-analytic study by 
iGranot et al.l (|2008l ). In this study 7-ray light curves pro- 
duced by an emitting thin shell in spherical expansion are 
computed at different energies taking into account the effect 
of the 77 opacity. Using our model and adopting the same 
prescriptions for the dynamics of the expanding shell and 
the properties of its emissio n (luminosity and s pectrum), 
we reproduce the results of iGranot et al.l (|2008l ) . The ex- 
cellent agreement between the two approaches is illustrated 
in Fig. O corresponding t o the example consid ered in the 
middle panel of Fig. 9 in IGranot et all l|2008h . The only 
difference occurs at very early times (first 0.1 ms) due to 
the limitation of the discretization used in our numerical 
method. The semi-analytical approach used bv lGranot et al.l 
(2008) does not suffer such a limitation but cannot be easily 
adapted to more complex outflows such as presented later 
in this paper . In th is example, the 77 opacity obtained by 
IGranot et al.l (|2008l ) and confirmed by our numerical calcu- 
lation is found to be a factor 92 smaller than the prediction 
from the single zone isotropic model, usi ng the normaliza- 
tion factor K^°^ from lAbdo et al.1 l|2009bl ) (see 3122]). This 
leads for instance to a cutoff energy in the time-integrated 
spectrum, defined as the energy -Ehe where r^-y ~ 1, which is 
~ 45r7ieC^ instead of 0.49mcC^. This strong discrepancy be- 
tween single zone models and detailed calculations confirms 
how crucial it is to take into account the precise geometry 
and variations of the radiation field when computing the 77 
opacity. 



(36) 3 THE 77 OPACITY IN INTERNAL SHOCKS 



In this section the model is applied to the dynamical evo- 
lution expected in the internal shock framework, where the 
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Figure 3. Opacity in a single p ulse: compari son with the 
semi-analytical work of Granot et alj (120080 • 77 absorbed 
lightcurves at four different energies e = EYiY^/m^c^ are plotted 
as a function of the observer time for a single pulse using the 
prescriptions and model parameters cor responding to the middle 
panel of Fig. 9 in iGranot et alJ KOOSh . Our numerical calcula- 
tion (solid line) is compared to the semi-analytical result (dashed 
line). Notations are the same (the observer time T and observed 
fluxes eFf/Fo have normalized values). The agreement is excellent 
except for T < 10~^, where the discrepancy is due to numerical 
resolution limitations (this corresponds to a true observer time 
tobs < 0.1 ms). 



This allows to follow the physical conditions in the shocked 
regions: Lorentz factor F,, conioving mass density p* and 
specific internal energy e« . A fraction e^ of the internal en- 
ergy is injected in relativistic electrons that are radiatively 
efficient. Therefore, the same amount of energy is radiated. 
We assume that each elementary collision produces a spher- 
ical flash as described in !j2]2]and that the comoving emis- 
sivity follows a broken power-law spectrum with low- and 
high-energy photon indices a and /3 and a peak energy 



E'p = KpUl 



(38) 



where K, x and y are constant. We use the values of x and 
y provided by the standard fast cooling synchrotron model 
with constant microphysics parameters, i.e. x = 1/2 and 
y = 5/2. In this case, the constant K depends on e^, on the 
fraction ea of the energy which is injected in the magnetic 
field and on th e fraction C of electrons that a r e acc elerated 
by the shock (iDaigne fc MochkovitchI Il998l . |2003| ). Once 
the dynamics and emission of internal shocks have been 
obtained, we apply the formalism described in Sect. [2] to 
compute the 77 opacity seen by high-energy photons. The 
assumption of a geometrically thin emitting region for each 
shock wave is well justifieqj, as electrons are in the fast 
cooling regime, limiting the size of the emitting region at 
the shock front to be of the order of c x t[^^, where t[^^ 
is the radiative timescale and is much shorter than the 
dynamical timescale idyn- We have tested that the results 
presented in this paper are independent of the choice for 
the time discretization of the dynamical calculation, as 
long as the time step is chosen to be short compared to 
the variability timescale Aivar in the outflow. The results 
described hereafter in the paper have been obtained using 
a time step tw/3000, except for the example shown in 
Fig. [8]where the time step has been decreased to tw/10000 
because of the short timescale variability in this particular 



whole prompt 7-ray emission is produced by electrons accel- 
erated by sho ck waves propagating within a relativistic vari- 
able outflow jRees fc Meszaroslll99i : iKobayashi et al.lll997l : 
iDaigne fc MochkovitchI Il998l ). A previou s study of the 7 7 
opacity in internal shocks was made by lAoi et al.l l|201(j ) . 
However the prescription used to compute r^-y was still ap- 
proximate, using the local physical conditions of the outflow 
where the high energy photon is emitted and applying them 
to an average formula of T.y-, as di scussed in H2.2.21 using the 
same normalization as in lLithwick fc Sari. (,2001, ). 



3.1 Description of the internal shock model 

We assume that the ejection by the central engine lasts for 
a duration iw and we consider the relativistic outflow at the 
end of the acceleration process, where it is assumed that the 
energy content is dominated by the kinetic energy. It is then 
described by the distribution of the Lorentz factor r(i) and 
of the kinetic power Ekin{t). Both quantities can vary on a 
timescale Atvar (variability timescale of the central engine), 
which leads to the formation of internal shocks at larger dis- 
tance. We model the dynamics of these internal shocks via 
a multiple shell model where the successive collisions be- 
tween shells mimic the propagation of shock waves withi n 
the relativistic outflow (|Daigne fc Mochkovitc"hlll998l . l200ol ). 



3.2 Use of the on axis 77 opacity as a good proxy 

From this point until the end of Sect. |3]the 77 opacity com- 
puted on-axis will be used as a proxy for the complete (lati- 
tude dependent) calculation. More precisely for a given flash, 
the approximate opacity 



T77 (Eue) — k 



(£HE;ee = 0) 



(39) 



will be applied to all high energy photons of energy i^HE 
radiated at the same radius, even those that are emitted 
off-axis. The constant factor k is determined below, from 
a detailed comparison with the exact calculation. This ap- 
proximation has the advantage of reducing the computing 
time without affecting the qualitative discussion. Even if 
the on-axis opacity (see Eq. (|18|) ) and the exact opacity r^-y 



•^ We have tested in the single pulse burst presented in Fig. 4 
that the contribution of a local isotropic radiation field associ- 
ated with the finite size of the emitting region is negligible as 
long as its width is smaller than r^ ct'^ ^/lOO. In the realistic 
case where this width is fixed by the synchrotron timescale, the 
additional contribution to the opacity is about five orders of mag- 
nitude below the main contribution. 
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Figure 4. 77 on-axis opacity as a good proxy. Same single 
pulse burst as in Fig. [7] (see Fig. [7] for the initial conditions in tiie 
outflow). Bottom panel: normalized piioton flux in the 260 keV 
- 5 MeV energy range as a function of observer time; Top panel: 
77 opacity seen by f GeV photons (source frame) as a function 
of the observer time. The solid line shows the flux-averaged value 
of the exact opacity (latitude dependent), as given by Ea. ll40ll : 
the short dashed line shows the opacity seen by photons emit- 
ted on-axis; the long dashed line shows the flux-averaged value 
of the approximate opacity defined by Ea. lj39|l using a constant 
factor fc = 3. Note that the curve for the on-axis opacity stops 
a-t iobs/(l -I- ^) — 1-33 s when the on-axis internal shock emission 
stops, whereas the flux averaged opacities are still well defined 
after this point, when the observed flux is only due to the high- 
(co-)latitude emission. The approximate value of Try-y (single zone 
model) is plotted as an horizontal thin dotted line using Ea. l|44|l . 
Note that this approximate formula has been calibrated using 
time-integrated spectra (see il3.3.2|l . 



(see Eq.(l7])) have different values, the ratio between the 
two remains fairly unchanged during the whole simulation 
as long as the on-axis emission is active. It is only at the 
very end of the prompt emission that the observed flux is 
dominated by the high-(co-)latitude emission and that this 
approximation is not appropriate anymore. 

To illustrate the validity of this approximation. Fig. [4] 
shows, in the case of a simple single pulse burslQ , the evolu- 
tion of the flux-averaged value of the exact opacity (latitude 



■^ In the examples presented in Figs.|4l|6l[71|8]and[9l high values 
of the kinetic power have been chosen as Fermj-LAT GRBs are 
among the brightest GRBs ever detected. For a jet opening angle 
from 1 to 10° , an isotropic equivalent kinetic power of 10^^ erg.s^^ 
corresponds to a true power of 1.5 • 10^^ to 1.5 ■ 10^"^ erg.s"^. 
In addition, the microphysics parameters are fixed to have an 
observed time-integrated peak energy -Ep,o = lMeV/(l -|- z) and 
the low- and high-ene rgy slopes are fixe d to the typical values 
a = —1 and /3 = —2.2 jPreece et al.ll2000ri in these five figures. 



X^i -^^(iobs, -Bhe) 



(40) 

where fobs is the observer time and the summation is done 
over all collisions contributing to the observed flux at iobs: 
for each i, Fi{toha, Ehe) is the contribution to the flux of 
collision i and T-yy{tobs,Qc,i) is the 77 opacity computed 
for collision i taking into account the emission colatitude 
Oc,i at time iobs- It is compared to the evolution of the 
opacity for on-axis photons and to the evolution of the flux- 
averaged value of the approximate opacity f.^^. deflned by 
Eq. (|39p using fc = 3. In this example, the 77 opacity evolves 
over more than 2 orders of magnitude whereas the ratio 
f-y^/T-y^ remains very close to unity, except at the end of 
the pulse, especially after the end of the on-axis emission. 
This small discrepancy is negligible in complex lightcurves 
as long as the central source is active. 



3.3 Minimum Lorentz factor in GRB outflows 

3.3.1 The case of GRB 080916C 

The first natural application of our model is the estimate 
of the minimum bulk Lorentz factor Fmin in GRB outflows, 
obtained from the constraint rT,^(i?HE,max) — 1, where 
i^HE.max IS the highest photon energy detected in the 
burst. To illustrate this aspect, we applied our approach 
to the case of one of the four brightest GRBs detected in 
the GeV range by Fermi, GRB 080916C. The results are 
shown in Fig. [S] Using our numerical model, a synthetic 
GRB is generated, which reproduces the main observa- 
tional features: the total radiated isotropic 7-ray energy 
(£iso = 8.8 X lO^"* ergs between 10 keV and 10 GeV), 
the spectral properties (i5p, a, j3 parameters of the Band 
functiorjfl), the envelop of the light curve and a short 
time-scale variability of 0.5 s in the observer frame. The 
study is focused on the most constraining time bin, bin 
'b', during which the highest observed photon energy was 
She, max = 3 GeV (16 GeV in the source rest frame): for 
this reason, only time bins 'a' and 'b' are reproduced in the 
synthetic GRB. These two intervals correspond to 32 % of 
the total radiated energy. The ejection by the central engine 
lasts for fw = 1.8 s, the injected kinetic power is taken 
constant and equals -Ekin(f) = 8 • 10^^ erg.s"'^ (isotropic 
equivalent value) and the assumed initial distribution of 
the Lorentz factor V{t) is plotted in the top-left panel of 
Fig. O This initial distribution leads to the formation of two 
shock waves. The flrst shock wave, propagating forwards, 
is short-lived as it reaches rapidly the front edge of the 
outflow. The second shock wave, propagating backwards, 
is the main source of radiation. Due to the initial shape 
of F(f), the physical conditions in the shocked region 
vary during the propagation, leading to two main pulses 
in the lightcurve, with additional variability on shorter 
timescales (see Fig. O top- right panel) . Despite its apparent 
complexity, this case remains relatively simple as there is 
only one propagating radiating front in the outflow. We 



^ We use the values given by lAbdo et al.l ll2009lj) for time bin 
'b'. 
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Figure 5. Minimum Lorentz factor for GRB 080916C. The two first panels are plotted for the critical case leading to T-y-|,(3 GeV) = 
1 in time bin 'b', i.e. for a mean Lorentz factor V = Fmin = 340. Upper left panel: initial Lorentz factor distribution in the outflow. Upper 
right panel: 7-ray lightcurvcs in the GBM/BGO band (8 keV — 5 MeV, top) and at 3 GeV (bottom). The lightcurves are plotted as a 



function of t„ 



-tn 



, where to^g trig is the observer time of the first detected photons. Lower left panel: evolution of t^^ at Ehe ■ 



GeV agains t the mean Lorentz factor in the outflow F, following our detailed modeling (solid line) and using the single zone isotropic 
model from lAbdo et al.l l|2009tJ) (dashed line). Lower right panel: time integrated spectrum over time bin 'b' for different mean Lorentz 
factors (the relative shape of the initial Lorentz factor distribution is kept the same) and reference spectrum without 77 annihilation 
(dashed line). 



find that the 77 opacity seen by high-energy photons is 
dominated by a local contribution, i.e. by the interaction 
with seed photons emitted almost at the same radius as 
the high-energy photon. Typically, we find that the opacity 
T77 (i?HE;Oo = 0) of a high-energy photon emitted on-axis 
at radius Re is built over a distance (. ~ (0.1 — > l)Rc- In 
addition, in this simple case, the geometry makes impossible 
head-on collisions: the interaction angle is necessary in the 
range < \1/ < l/F, resulting in a seed radiation field which 



is highly anisotropic. Therefore, the situation considered to 
model time bins 'a' and 'b' in GRB 080916C is not very far 
from the simplest case of a spherical flash, and, as discussed 
below, the approximate formulae derived in H2. 2. 21 remains 
relevant even in this more complex case. 



The minimum mean Lorentz factor Fmin is obtained by 
requiring that r.y^ (i^HE.max) < 1 (see Fig. [S] lower panel). 
With the detailed calculation, we find a minimum mean 
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Lorentz factor Fmin = 340, i.e. a factor 2.6 lower than the 
value Fmin = 887, which wa s obtained from an approximate 
single zone isotropic model iJAbdo et al.ll2009bf l. Even more 
remarkable, the whole initial distribution of the Lorentz 
factor used in this model of GRB 080916C (from 170 to 
700) remains below the "minimum" value of the Lorentz 
factor derived from single zone models (see Fig. [5] upper 
left panel). 

Compared to GRB 80916C, the case of GRB 090926A 
IjAckermann et al.l l201ll ) is more complex : the time- 
integrated high-energy spectrum shows a cutoff at about 
1.4 GeV. If this cutoff is interpreted as due to 77 anni- 
hilation, it allows a direct measurement of the Lorentz 
factor instead of a minimum value. It is however difficult 
to confirm this interpretation from the spectral shape of 
the cutoff and other physical explanations are possible 
(see H3.4|l . Another source of complexity is due to the 
presence of an additional component at high-energy, which 
is not necessarily produced in the same region than the 
soft gamma-ray component (see the discussion in Sect. H}. 
Nevertheless, interpreting the cutoff a s a signature of 77 
annihilation, lAckermann et al.l l|2011f) obtain a Lorentz 
factor of F ~ 720 ± 76 with the single zone model, and 
a value F « 220 (i.e. a factor ~ 3 lower) with a more 
realistic approach b ased on the formalism developed by 
iGranot et all (|2008l ) (see ig^] and tj3.3.2p . Using the ap- 
proximate formula given in Eq. (|54[) below, which has been 
calibrated on the numerical results presented in this paper 
(see ij3.7.3|l . we obtain a very similar value F ~ 250. 

If the compactness argument imposes a minimum 
Lorentz factor for the outflow, we note on the other hand 
that - if the prompt emission has an internal origin - there 
is also an upper limit on the Lorentz factor, imposed by the 
condition that the internal dissipation should occur before 
the deceleration by the external medium becomes signifl- 
cant, i.e. at radii lower than the deceleration radius: 



Ra 



E-kin 



47r 



F ylc2 



(41) 



where F is the mean Lorentz factor and _Ekin the total kinetic 
energy of the outflow, and A and s parametrize the density 
profile of the external medium which follows 



A 
Pcxt - -^ ■ 



(42) 



condition _Ris < -Rdoc imposes F < 420 in the wind-like 
medium (resp. F < 540 in the uniform medium). Assum- 
ing dense external media such as expected around a massive 
star, previous estimates of Fmin ~ 887 were not consistent 
with this constraint, whereas the more accurate estimate of 
Fmin = 340 obtained in the present study fulfills the condi- 
tion for both types of external medium. 

3.3.2 Minimum Lorentz factor: detailed modeling vs 
simple estimates 

The difference in the estimate of Fmin from single zone 
isotropic models and the detailed approach presented here 
has two origins: as shown in H2.2.2I the normalization 
obtained with single zone isotropic models is larger than 
the normalization obtained by considering the interaction 
of a high-energy photon with the anisotropic field created 
by a single flash, leading to a ratio in the estimate of the 
minimum Lorentz factor equal to [Kg'^^ /Kq) ^Ci-") ~ I.5 

llAbdo et al.ll2009bh (r esp. {K^^°^ /Ko)^^^ ~ 1.3 ^ 1.5, 
iLithwick fc Saril 120011 ) for /? = -2 ^ -3. However the 
difference illustrated by the example of GRB 080916C 
is much larger, due to the fact that the dynamics of the 
outflow is properly taken into account in the detailed 
approach presented here. Such effects are important, see 
e.g. the dependency on Rb in Eq. (|25p . 

Note that a precise comparison between single zone 
models and our detailed approach is difficult, as the value 
of the Lorentz factor is not unique in our case. Indeed, the 
variations of F are at the origin of the prompt emission, 
as we use the internal shock framework. In the example of 
GRB 080916C, the value Fmin = 340 is the mean value for 
a distribution varying from Finf = 170 to F^up — 700. In 
the approximate formula for r-f-, obtained in ^2.2.21 the de- 



pendency T~f.y OC Fq 



2(1-/5) 



can be seen as the product of 



Fq coming from the dependency of Tj^, on the emission 
radius, and Fq (i.e. F^^* for /3 = —2.2) coming from 

the dependency of r^.-, on the photon field. In the inter- 
nal shock model, the radius is mainly fixed by the mini- 
mum Lorentz factor Finf within the outflow. Indeed, in a 
two shell collision with a ratio k — Fsup/Finf between the 
Lorentz factors and a variability timescale Atvar, the inter- 
nal shock radius depends only on Finf for efficient collisions. 
It equals R ~ 2Ffnf cAtvar with a precision better than 10 % 
for K > 3.2. On the other hand the Lorentz factor relevant 
for the radiation field should be the Lorentz factor in the 
shocked regions, which is not too different from the mean 
Lorentz factor F in the outfiow. This leads to 



'77 "^ ^ inf -■- 



As Rdec decreases as F ^""^ , this condition may become 

severe for high values of Fmin. This is especially true for 

the internal shock model, where the typical dissipation ra- 

2 

dius increases with the Lorentz factor as F . In the case 

of the synthetic burst used to model bins 'a' and 'b' of 
GRB 080916C, a typical efficiency of 5 % for the inter- 
nal shock phase leads to a total kinetic energy of the 
outflow -Ekin — 1.4 ■ 10^® erg. Therefore, for a wind- like 
(resp. uniform) medium with A = 5 ■ 10^^g-cm~^ and 
8 = 2 (resp. A = lO^nip cm~^ and s = 0) , the deceler- 
ation radius equals -Rdcc ~ 1.9 • 10^^ (F/lOO) cm (resp. 

Tidec = 5.6-10" (F/lOO) "^^^ cm). On the other hand, inter- 
nal shocks are found in the simulation to occur over a radial 
range Ris ~ (8.7 ■ 10^^ -> 6.1 ■ lO") (F/lOO)^ cm. Then the with K'o = CiA'o = Ci {2^+'^'^l{l3)) 



(43) 



All the examples presented in this paper agree well with 
this scaling law. We flnd that a good approximate formula 
for the 77 opacity can be deduced from the comparison of 
Ea. (|25|) and the detailed numerical calculation presented in 
this paper and is given by 



i{Ee 



K 



AootE-c, 



Atv (cAf 



var J -C/| 



P,0 \ 



/ (moC^) 



2\2\ 1+/3 



-Bhe-Bp,o 



4f^2(l+;9) 
•- inf J- 



(44) 
The numerical 



Do Fermi-LAT observations imply very large Lorentz factors in GRB outflows ? 11 



constant Ci ~ 4 • 10~^ has been obtained from the com- 
parison of Eq. (|44p with the result of the fuU numerical 
calculation in a series of runs based on the synthetic single 
pulse GRB shown in Fig. [4] In this example, we have 
Tint = 100, r = 290, -Bp.o = IMeV, /? = -2.2, Aivar ^ 0.4s 
and the radiated energy above _Ep,o is 1-3 ■ 10^'^ erg, i.e. 
frad = 5.2 • 10*^ erg after a correction by a factor 0.4s/ls 
to estimate f'rad over Afvar. Then each of these parameters 
has been varied in different series of runs, to allow for the 
calibration of Ci . As seen in Fig. [l] it is not relevant to com- 
pare the approximate value of the opacity given by Eq. 1)44^ 
with the result of the time-dependent calculation, where 
the opacity varies by several orders of magnitude during the 
evolution of the pulse. Therefore, the comparison is based 
on the time-integrated spectrum over the full duration of 
the burst. A "time-averaged" opacity is computed from 
the ratio of the unabsorbed over the absorbed spectrum, 
and then compared to Ea. H44p . As the single zone model 
predicts an exponential cutoff at high-energy, whereas 
the detailed calculation shows a more complex shape (see 
ij3.4p . the comparison can be done only in the relevant 
energy range, i.e. in the region where the opacity is close 
to unity. At higher energy, the 77 opacity is overestimated 
by the single zone model. At lower energy, we find that the 
discrepancy on the 77 opacity is less than a factor of 2, 
which would lead to a factor 2^'^'^^^~^'' for the minimum 
Lorentz factor, i.e. about 10 %. In the series of runs which 
were used to calibrate Ci, we have varied the mean Lorentz 
factor r from 100 to 1000, the ratio k = Fsup/Finf from 2 
to 8, the kinetic power Skin from 10^° to lO^^erg.s"^, the 
duration tw from 0.1 to 10 s, and the peak energy _Bp,o 
from 100 keV to 10 MeV and we have found that - after 
calibration - the discrepancy for the 77 opacity remains less 
than a factor of 2, except for low k, where it can increase 
to a factor of 10. This is because the simple scaling law 
between the internal shock radius and the Lorentz factor 
Finf is not valid anymore for k < 3. We have checked a 
posteriori that the discrepancy for the opacity is less than 
a factor of 2 for all the cases presented in this paper. 



T o compare to the esti mates by lAbdo et al.l (J2009al l 
(resp. iLithwick fc Sarill200ll ) where the details of the dis- 
tribution of the Lorentz factor are unknown, we use Eq. (|44[) 
with Finf = F. This leads to a minimum Lorentz factor 
Fmin which is reduced by a factor {K^°^ / Kq)"^^^ ~ 2.4 

(resp. {Ky^°^ / Kq)"^^^) ~ 2.2) for /3 = -2.3. This is in rea- 
sonable agreement with the pre-Fermi study presented by 
iGranot et al.l (|200g ). In this earlier work, a semi-analytical 
calculation of r-y^ using a simple prescription for the propa- 
gation of the shock wave responsible for the pulse (constant 
Lorentz factor) but taking into account the exact radiation 
field led to a normal ization in the equiv alent of Eq. (|44|) given 
by (see Eq.(126) in lCranot et al.ll2008l '): 



^/G08 
^0 



1.86 • 10" 



C2 



'P 



-5/3 



(45) 



where C2 — 1. For /3 



-(3) 



factor [k^^'/K'''') '''-' 
minimum Lorentz factor 



-2.3, this leads to a reduction 
2.6 for the estimate of the 



3.4 Shape of the 77 attenuated time-integrated 
spectrum 



As pointed out bv lGranot et al.l l|2008l ). due to the temporal 
evolution of t^^, the opacity cutoff in a time-integrated 
spectrum will be smoother than a sharp exponential decay: 
the cutoff transition will be close to a power-law steepening. 
The detailed model presented here is appropriate to char- 
acterize from an observer point of view this time evolution 
effect. The smoother r^-y transition is due to the fact that 
the time integrated spectrum is a superposition of instant 
spectra which can have different 77 cutoff photon energies. 
This time evolution takes place within a given 7-ray pulse, 
and can be even stronger in a complex burst w here the light 
curve is made of many pulses IjAoi et al.ll20ld ). In addition, 
in the latter case, cross-interactions between different pulses 
can become important and strongly influence the time 
evolution of the opacity: such cross-interactions are fully 
taken into account for the first time by our approach. 



To illustrate this effect, two examples of synthetic 
GRBs are shown in Fig. [G] In the case of the "mono-pulse" 
GRB the 77 opacity seen by high energy photons evolves 
smoothly, with a regular decay (except at the very begin- 
ning when t he interacting phot on field progressively builds 
up, see also iGranot et al.l 120081 ): as a result, even for time 
intervals corresponding to a noticeable fraction of the burst 
duration, the 77 cutoff in the time-integrated spectrum 
remains close to an exponential cutoff. In the case of the 
"complex" multiple-pulse GRB, the 77 opacity shows much 
stronger variations, due to the intrinsic variability of the 
outflow, and consequently, the 77 cutoff is much closer to a 
power-law steepening. 



To our knowledge, GRB 090926A is the best case where 
a cutoff in the spectrum a t high energy has been cl early iden- 
tified at about 1.4 GeV IjAckermann et al.ll2011h . However 
the shape of this cutoff is weakly constrained due to the poor 
photon statistics at high energy. As a result, it is difficult 
to further discuss if this break is due to the 77 annihila- 
tion or some other spectral effect. For example in a leptonic 
scenario, the main Band component could be produced by 
synchrotron emission whereas a second sub-dominant com- 
ponent at higher energy could be due to in verse Compton 
scatterings in Klein Nishina regime (see e.g. iBosniak et al.l 
[2OO9I). In this scenario, it is difficult to distinguish the cut- 
off due to 77 annihilation from the expected cutoff of the 
inverse Compton component at high energy, without a good 
characterization of the spectral shape above 1 GeV. Even in 
the absence of clear spectral signatures of a cutoff for most 
Fermi GRBs, an evidence for an attenuation at high-energy 
in a large fraction of Fermi bursts is given by the compar- 
ison of the LAT detection rate and the theoretical rate ob- 
tained by extrapolating the G RB spectrum measured by the 
GBM (|Beniamini et al.ll201H ). Unfortunately, this method 
does not allow any comparison of the spectral shape of the 
cutoff with theoretical predictions. 
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Figure 6. Temporal and spectral evolution of the 77 opacity - Example of two synthetic GRBs. The upper part of the four 
figures illustrates the case of a mono-pulse GRB with no additional temporal substructures whereas the bottom part show the case of 
a complex multiple-pulse GRB. The Lorentz factor during the ejection rises monotonously from 100 to 400 in the first case, leading to 
a collision between the slow and the fast parts of the outflow, and from 100 to 700 in the second case with the addition of variability 
on a shorter time-scales. The kinetic power of the outflow is assumed to be constant, with -Ejjin = lO^'^erg.s"^ in the first case and 
5 ■ 10^^ erg.s"^ in the second case (so that the radiated energy 2 • lO^'^ ergs is the same in both cases). Upper left panel: initial Lorentz 
factor distribution of the outflow as a function of the injection time ti^j. Upper right panel: the evolution of the approximate opacity 



T-y-|, , as defined in ijc{.2l seen by 1 GeV photons (source frame) as a function of the observer time t^ 



''obs, trig 



The curve is plotted only 



as long as the on-axis emission is active. Lower left panel: the 7-ray light curve at 1 MeV (source frame) as a function of the observer 
time toi,g — tobs,trig) where tobs,trig is the observer time of the first detected photons in the complex GRB case. Three integration time 
bins for the simulated spectra (see lower right panel) are delimited by vertical dotted lines and labeled 'a', 'b' and 'c'. Lower right 
panel: Time-integrated spectra corresponding to time bins 'a', 'b' and 'c' (dotted line) and 'a-fb~fc' (dashed line). In the latter case, the 
time-integrated spectrum without 77 absorption is plotted in solid line. 



3.5 Is the delayed onset of the GeV emission a 
signature of the 77 opacity ? 

3.5.1 Effect on the htgh-energy lightcurve of a time 
evolving 77 opacity 



The high energy emission (above 100 MeV) detected by 
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Figure 7. Delayed onset produced by the 77 opacity. We consider a single pulse burst generated by internal shocks in a relativistic 
outflow. Upper left panel: initial Lorentz factor distribution of the outflow as a function of the injection time iinj . The kinetic power is 
constant and equals lO^^erg.s"^. Upper right panel: the approximate opacity t^-/, as defined in i|3.2l seen by 1 GeV photons (source 
frame) is plotted as a function of the observer time tobs- Lower right panel: high-energy light curves calculated for She = 10, 100, 300 
MeV, 1,3, 10, 30 GeV (source frame). The delayed onset begins to show up above 1 GeV, which would correspond to 300 MeV for 
a GRB at z = 2 (such as GRB 090926A) or 200 MeV for a GRB at 2 = 4 (such as GRB 080916C). Lower left panel: time-integrated 
spectrum over the whole duration of the burst. The unabsorbed spectrum is plotted in dashed line. The vertical dotted lines mark the 
photon energies at which light curves are calculated. 



compared to the softer 7-ray emis sion (below 5 MeV). 
The analysis by IZhang et al.l (|201l|) indicates that such 
a delayed onset is present in at least 7 in a sample of 17 
GRBs detected by Fermi-hAT . This feature seems to be 
common to long and short GRBs and its origin is debated 



(jGranot et al 



2OI0I). Am on g the proposed explana tions 



(see e.g. IZou et al.M2009l : 111 120101 : iToma et all |2009| ). the 
possibility that this delayed onset is induced by a 77 opacity 
temporal evolution effect has already been discussed by 



lAbdo et al.l l|2009bl ): as the shock wave producing the 7-ray 
emission expands to larger radii, the opacity seen by the 
high energy photons evolve from an optically thick to an 
optically thin regime. The model developed in the present 
study is well appropriate to investigate this possibility 
in more details. Fig. [7] illustrates a simple example of a 
synthetic burst displaying a delay between the high-energy 
and soft 7-ray lightcurves due to 77 annihilation. The 
initial Lorentz factor distribution of the relativistic outflow 
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is simple, with the formation of two shocks when the slow 
part of the ejecta catches up with the fast part: an internal 
"forward" shock sweeping the slow part and an internal "re- 
verse" shock sweeping the fast part. The internal "forward" 
shock quickly disappears and contributes to the observed 
7-ray emission only at early times (iobs/(l + z) < 0.21 
s). On the other hand the internal "reverse" shock has a 
longer duration; it forms at i? = 8 ■ 10^'^ cm and propagates 
until 7? = 3 • 10^^ cm, leading to most of the prompt 7-ray 
emission. During the propagation of the shock, two effects 
contribute to a reduction of the 77 opacity for high-energy 
photons detected at later times: both the radius and the 
Lorentz factor of the shocked material are increasing. 

This example shows that the temporal evolution of the 
77 opacity can actually induce a significant delay between 
the high energy and the soft 7-ray emission. The synthetic 
burst used in Fig. [S] to model bins 'a' and 'b' of GRB 
080916C gives another example of a delayed onset at 3 GeV 
induced by an evolving 77 opacity. The first pulse is pro- 
duced at lower radii and in lower Lorentz factor material 
and is therefore strongly absorbed. For this reason, it is al- 
most suppressed in the 3 GeV lightcurve, whereas the second 
pulse is well visible. 



1 GeV does not allow to strongly constrain this potential 
behavior. 

In a simplified picture one could also expect that the 
delay could not be larger that the duration of the 7-ray 
pulse. Within the internal shock model, this would be 
true if each 7-ray pulse was due to a distinct shock wave. 
This is not necessarily the case as the observed lightcurves 
often show many superimposed pulses, which suggests that 
another, more complex, situation is possible, where a single 
propagating shock wave is responsible both for the envelope 
of the lightcurve and the superimposed short duration 
pulses, due to modulations in the kinetic power and/or the 
Lorentz factor of the outflow. For instance, the synthetic 
GRB shown in Fig. [5] to model bins 'a' and 'b' of GRB 
080916C exhibits two main pulses with a superimposed 
short time-scale variability and is nevertheless due to a 
single shock wave with variations in the shocked region due 
to the modulation in the Lorentz factor (see top left panel 
of Fig. [Sj. The delayed onset in the GeV lightcurve induced 
by the evolving 77 opacity is of the order ~ 5 s, i.e. much 
larger that the typical duration of short duration pulses in 
this burst (~ 0.5 s). Such a delay is in good agreement with 
_Fermi-LAT observations. 



3.5.2 Characteristics of the 77 induced delay 

Consistency with the spectral analysis. Testing this 
theoretical idea for the observed delay remains difficult 
as the precise specifications of the high-energy 7-ray 
spectrum in Fermj-LAT GRBs are not fully understood. 
In some cases, the observed time integrated spectra are 
consistent with a unique Band function component over 
the GBM and LAT energy ranges (see e.g. GRB 080916G, 
lAbdo et al.ll2009a ). but in some other cases, they show an 
additional component at high energy ( e.g. GRB 090510, 
lAckermann et all l20ld GRB 090902B lAbdo et all l2009al 



and GRB 090926A. lAckermann et aljboill . It is then 
unclear whether the observed evolution of the high-energy 
slope P, such as the steepening from time bin 'a' to time bin 
'b' in GRB 080916C, is real or related to the evolution of 
an additional high-energy weak component. In the former 
case, this spectral evolution could also contribute to the 
observed delay. Note that even if /3 is intrinsically constant, 
the time evolution of the 77 opacity discussed here to 
explain the delayed onset of the GeV emission would also 
affect the measured value of the high-energy slope, due 
to the power-law steepening of the spectrum (see tj3.4p . 
Clearly, despite a very significant improvement compared 
to previous instruments, Fermi-hNT cannot yet provide a 
detailed picture of the high energy part of the spectrum 
where several effects are expected to superimpose: shape 
of the high-energy tail of the main component, location 
and shape of the additional weak high energy component, 
location and shape of the 77 attenuation. 

Predicted duration of the delay. An interesting pre- 
diction of the 77 induced delayed onset of the .Ffermi-LAT 
emission is that the delay should increase when observing 
at higher energy, since the transparency regime is reached 
at larger radii (and so at larger observer times) for photons 
of higher energies. However the paucity of photons above 



A high-ene r gy p r ecursor ? It has been argued by 
iGranot et all (J200a ): lAbdo et all (|2009bl ) that a high en- 
ergy emission should also be detected at the very beginning 
of the burst while the photon field inducing the 77 opacity 
is still building up. However this early transparency phase is 
very short (of the order of At obs — 100 ms for the synthetic 
burst presented here) and occurs when the 7-ray emission 
is still at the beginning of its rise. As seen in the bottom 
right panel of Fig. (T) the high-energy precursor is visible at 
1 GeV (source frame) but remains extremely weak. The rela- 
tive amplitude of the precursor and the main pulse increases 
at higher energy but the flux becomes very low. It is then 
unlikely that such a high-energy precursor can be detected 
by current instruments. 



3.6 Temporal smoothing at high energy 

Within the internal shock model, another expected signature 
of the 77 annihilation process is the smoothing of the short 
time-scale variability in the high energy light curves. Indeed, 
in complex bursts where the 7-ray emission is made of sev- 
eral shock waves, the opacity at a given frequency will be 
larger for photons emitted at smaller radii where the small- 
est time-scales of the observed light curves are produced. 
This effect is illustrated in Fig. [8] where we have simulated 
a complex multi-pulses GRB, using a highly variable dis- 
tribution of the initial Lorentz factor in the outflow. The 
variability has been produced on different timescales from 
100 ms to 4 s, with a total duration of the ejection of 30 s. 
Lightcurves are plotted at three energies 10 MeV, 300 MeV 
and 3 GeV (source frame). The attenuation due to 77 anni- 
hilation becomes strong above a cutoff energy which evolves 
strongly during the burst. In particular, short pulses tend 
to have a lower cutoff energy, because they are produced 
at smaller radii. The cutoff remains however always above 
10 MeV, so that the first lightcurve is un-attenuated and 
shows all the variability initially introduced in the Lorentz 
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Figure 8. Temporal smoothing effect of the 77 opacity at high energy. A complex multi-pulse burst has been generated by 
assuming an outflow ejected for 30 s with a constant kinetic power 2 • 10^^ erg.s^^ and an initial Lorentz factor varying between 50 and 
350 with modulations on timescales varying from 100 ms to 4 s. Left: initial distribution of the Lorentz factor in the outflow. Center: 
time-integrated spectrum. The unabsorbed spectrum is plotted in dashed line. The three vertical dotted lines indicate the energies 
corresponding to the three lightcurves in the next panel. Right: lightcurves at 10 MeV (top), 300 MeV (middle) and 3 GeV (bottom) 
(source frame). Note the attenuation of the short duration pulses in the 3 GeV lightcurve. 



factor distribution. At higher energies, the lightcurves show 
less and less variability, the suppression affecting mainly the 
short duration pulses. It has been proposed that the study 
of the variability in the GeV lightcurve could be used to test 
whether GeV photons observed by .Fermi-LAT have an in- 
ternal ori gin, or are produced by th e exte rn al shock as pro- 
posed by iKumar fc Barniol DuranI (|2010l ) ; iGhisellini et al.l 
l|201Cl) . In practice, the temporal smoothing effect described 
here could make such a test difficult. 



3.7 Other sources of opacities in the outflow 

Thomson scatterings by leptons in the outflow are an ad- 
ditional source of opacity to take into account in GRBs. 
Both primary electrons injected in the outflow by the cen- 
tral source with baryons, and secondary electron-positron 
pai rs produced by 77 ann ihilation should be considered (see 
e.g lLithwick fc Sari|[200ll ). 



3. 7. 1 Thomson opacity from primary electrons 

The Thomson opacity Te due to the primary electrons can 
be compute d in a variable outflow followin g the procedure 
described in lDaigne fc MochkovitchI (|2002h to estimate the 
photospheric emission expected in the internal shock model. 
Fig. [9] presents the result of such a calculation for a syn- 
thetic GRB very close to the example used in Fig. [T] except 
for a slightly higher Lorentz factor (mean value of 360 in- 
stead of 290 chosen to have well separated contributions to 
the opacity in Fig. [9|. The 77 opacity at 1 GeV (source 
frame) is plotted as a function of the observer time and fol- 
lows the evolution which has already been described earlier. 
The Thomson opacity due to primary electrons continously 
decreases as a function of the observer time, due to the di- 
lution associated to the radial expansion of the radiating 
shocked region. As expected, if we exclude the very begin- 
ning of the internal shock phase, most of the internal shock 
propagation occurs well above the photosphere of the out- 
flow (tc <^ 1). Note that this calculation is not as detailed 



as the calculation of the 77 opacity presented here, as we 
limit it to on-axis photons ( for a detail ed discussion of the 
off-ajds photosphere, see e.g. |Pe'eiil2008l ). Another source of 
uncertainty is the initial number Ye of electrons per nucleons 
in the outflow. It has however a limited impact, at most a 
factor 2 in Te. Everywhere in this paper, we adopt Ye = 0.5. 

3.7.2 Thomson opacity from pair-produced e 

The detailed calculation of the opacity t± due to secondary 
pairs of leptons produced by 77 annihilation is beyond the 
scope of this paper. We estimate r± with several simplify- 
ing assumptions : (i) we assume that the Thomson regime 
is valiqj which is reasonable since even if the pairs e could 
be highly energetic right after their production, they will ra- 
diate rapidly their energy. In practice, pairs will be mainly 
produced by the annihilation of photons at energy Scut de- 
fined by T~f.y (Ecut) — 1, with seed photons interacting with 
a typical interaction angle ^ ~ l/F at the threshold of pair 
production, i.e. at energy 

(2Fmec2)' 



Ecu 



Ec- 



(46) 



The two pair-produced leptons in this case have an energy 
(lab frame) 



T^±mcc 



Ecut + Ecut 



(47) 



For Ecut ^ 300 MeV to 1 GeV, and F ~ 300, this leads to low 
values of the Lorentz factor of the leptons, 7± ~ 2 — 4. The 
scattering of typical 7-ray photons at energy Ep^ (source 
frame) by these leptons will occur in the Thomson regime if 
the photon energy in the proper frame of the lepton is small 
compared to the rest mass energy of the lepton, i.e. if 

7±-Ep,o 



w± = 



FtTIoC^ 



<1. 



(48) 



° If the Klein Nishina corrections were not negligible, the inter- 
action cross-section would be smaller, reducing t±. 
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Figure 9. Comparison between the 77 and the Thomson opacities. We consider a single pulse burst generated by an outflow 
with a distribution of Lorentz factor from 125 to 500 with a similar shape as in Fig. [T] a total duration of 1 s and a kinetic power of 
^kin = 10^^ erg.s"^. Left •panel: opacities seen by photons emitted on the line of sight as a function of the observer time 1^^,^'- 77 opacity 
T^^ seen by 1 GeV photons (source frame) in solid line; Thomson opacity To due to the primary electrons in dashed line; Thomson 
opacity t-j- due to pair-produced e^ in dotted line. These three opacities are plotted only as long as the on-axis emission is active. Right 
panel: the cutoff energy due to 77 (source frame, top) and the Thomson opacity due to primary electrons (dashed line, bottom) and 
secondary leptons (dotted line, bottom) are plotted as a function of the mean Lorentz factor F, keeping the same relative shape for the 
initial distribution of the Lorentz factor in the outflow. The evolutions of Ecut and t± show a break at F ~ 170 when Ecut = itself (see 
Eg. HSOI l). In addition, the evolution of Ecut is interrupted at F ~ 810 when Ecut = ^p,0- The location of the case shown in the left panel 
with F = 360 is indicated with a vertical thin solid line. The minimum value of the Lorentz factor for the outflow to remain optically 
thin for the Thomson opacity (here dominated by the contribution of secondary leptons) is indicated with a big dot and corresponds to 
F = 270. 



This condition is clearly fulfilled for the typical values of y± 
estimated above; (ii) it is assumed that pairs are produced 
at the location where the photons eventually annihilating 
were emitted. This approximation is excellent for simple 
GRBs where 77 cross-interactions between different shock 
waves are limited, since in such simple cases most of the high 
energy photons annihilate very close to their emission loca- 
tion with seed photons produced by the same shock wave. In 
more complex situations, this approximation is less justified. 



The value of t± is plotted as a function of the observer 
time in Fig. [5] As expected, it decreases with time, due both 
to the decrease of the 77 opacity leading to a decrease of the 
pair production, and to dilution due to the radial expansion. 



Note that we do not include in this study the contri- 
bution to radiation of secondary leptons, such as pairs pro- 
duced by 77 annihilation, or additional pairs produced by 
cascades if the first generation of pairs produce new high en- 
ergy photons that can annihilate (the cascade process would 
stop when the secondary photons have too low energies to 
produce new pairs). As the fraction of the radiated energy 
in annihilating photons is low in all examples presented here 
(the 77 cutoff being at high energy), the neglected contri- 
bution remains weak. 



3.7.3 What is the dominant opacity ? 

In the example shown in Fig. [51 all opacities (77 opacity at 
1 GeV and Thomson opacities due to primary electrons and 
pair-produced leptons) are decreasing with time. For most 
of the burst, especially at the maximum of the lightcurve 
at tobs — 0.4 s (see Fig. [7}, the dominant opacity is the 77 
opacity at 1 GeV, which means that, assuming that a GeV 
photons is detected in this burst, the most constraining 
limit on the minimum Lorentz factor in the outffow will be 
provided by the 77 opacity limit discussed earlier. However, 
the dependency of these three opacities with the Lorentz 
factor of the outfiow is not the same. Then, depending on 
the maximum energy of detected photons, the limit could 
be provided by any of the three constraints, as discussed in 
iLithwick fc Saril l|200ll ). 

In the right panel of Fig. |§] we have plotted the value 
of the Thomson opacity due to primary electrons and sec- 
ondary leptons and the cutoff energy due to the 77 opac- 
ity as a function of the mean Lorentz factor in the outflow 
(other parameters are kept constant: see caption). We find 
that these values are in good agreement with the following 
approximate formulae: 

(i) Opacity due to 77 annihilation: Eq. (|44p : 
(ii) Thomson opacit y due to primary electrons: from 
iDaigne fc Mochkovitchl l|2002l ). this opacity is approxima- 
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tively given by 



C2 






p-2p-3 

-■- int-"- : 



(49) 



where C2 — 0.2 is a correction factor obtained by comparing 
tiie approximate formula and the numerical calculation, and 
f-y is the efficiency of the prompt emission phase, so that 
the kinetic energy of the flow is given by Ermd/f-y- For the 
example used in Fig. [Ul this efficiency equals /^ ~ 0.02; 

(iii) Thomson opacity due to pair-produced leptons: 
this opacity can be estimated in a similar way as in 
iLithwick fc Saril (|200fh . taking into account the correction 
factors for the 77 opacity: we define E^^it = 2rmcC'^ as the 
energy (source frame) of photons that self-annihilate with 
an interaction angle ^ = f/F. Then the Thomson opacity 
due to pairs is approximatively given by 



T± 



with 



C3X 



l+,8 
2 



-.r-F^ 



1+" p 



2p/3 



T« 



AoCTTSn 



47r (cAtvar) £| 



p,0 



if -Ecut > -Escif 
if -Ecut < Escif 

1 + /3 



E, 



P,0 



(50) 



(51) 



and where C3 ~ 3 is a correction factor introduced from the 
comparison with the detailed calculation. The cutoff energy 
due to 77 annihilation can be estimated from Eq. (|44p : 



-Eci 



m-oC (KoTt) 



pTT?p^ 
^ inf ^ 



The condition Ecut > Escif is equivalent to 



K 

21+/3 



4-pl+/3 



(52) 



(53) 



The approximate formulae for r-y^, and t± assume a power- 
law spectrum, as in H2.2.2I Therefore, these equations are 
valid if -Ep,o is understood as the energy above which the 
power-law with slope j3 is observed, frad is the radiated 
energy above -Ep,o over a timescale Atvar and E^cut and 
Sseif are above -Ep,o. For a comparison with Fermi obser- 
vations, frad can be estimated from the measured photon 
fluence above _Ep,o (see Eq. (|26p ). The variability timescale 
can be estimated from the observed variability timescale 
in the GBM lightcurve. In addition, corrections due to 
the redshift z of the source should also be taken into account. 

We have checked that all simulations presented in this 
paper (including the model of GRB 080916C in Fig. [5} cor- 
responds to situations where the GRB outflow during the 
emission is optically thin for Thomson scatterings by pri- 
mary electrons or pair-produced leptons. However, from the 
scaling laws listed above, it appears that the limit on the 
minimum Lorentz factor in the outflow could be constrained 
by any of the three opacities, depending on the highest 
energy photon detected in a given burst (jLithwick fc Saril 
[2OOII). With our corrected estimates of the three opacities, 
we obtain the following limits, assuming F ~ Finf ~ F to 
allow a comparison with previous single-zone studies: 

(i) Opacity due to 77 annihilation. From Ea. (|44|) . if the 
highest energy photon has an energy Emax (source frame). 



the minimum Lorentz factor is given by 

yloO"Tfrad 



r ■ 

-L mm,- 



CiKo 



ATv{cAt^^r) Ep,0 



„2\2 \ 2 1-/J 



-L/max-L/p,0 



(54) 

If the cutoff energy _Ecut is identified, the previous equa- 
tion can be used with -Ecut instead of -Emax to estimate the 
Lorentz factor; 

(ii) Thomson opacity due to primary electrons. From 
Eq. (|49|l . the minimum Lorentz factor for transparency is 
given by 



r ■ 

-*- min,( 



C2 



YcQ-tE^-s 



SirrUpC'^ (cAfvar) A 



(55) 



(iii) Thomson opacity due to pair-produced leptons. From 
Eq. (|50p . the minimum Lorentz factor for transparency is 
given by 



r ■ 

-L mi 



n,± 



<-3Ci^j r, if^-R^o< — TW 

/„^l + ^\3^ 3^ cCi jv- . 2^ + ^-^ 



(56) 

We have plotted in Fig. [TO] the three constraints in a 
plane Lorentz factor - radiated energy, based on Ea. (|54|) . 
Eg. ((55)) and Eq.(l56|. In panel (a), for /3 = -2.2 and -Ep,o = 
IMeV (source frame), the minimum value of the Lorentz 
factor is determined by the 77 opacity as long as the cutoff 
energy is above ~l-3 GeV (the precise limit depends on the 
radiated energy frad). For a reasonable range in the radi- 
ated energy, f rad = 10*^ — > 10^^ erg, and for cutoff energies 
below 100 GeV, this leads to minimum Lorentz factors from 
21 to 440. If the cutoff energy due to 77 annihilation is be- 
low 1 GeV, then the minimum Lorentz factor is determined 
by the Thomson opacity due to leptons. In the present case, 
it is dominated by the opacity from pair-produced leptons. 
This limit goes from F = 18 for f^ad = lO"^ erg to 250 for 
f rad = 10^^ erg. Notice that the cutoff energies listed here 
are measured in the source frame. In the observer frame, the 
limit at 1-3 GeV is located in the low-energy part of the LAT 
spectral range (100 MeV-1 GeV) for usual GRB redshifts. 
In panel (b), we study the impact of the high-energy spectral 
slope (3 by using /3 = —2.5 and -Ep,o = 1 MeV (source frame). 
The main effect is to increase the region where the 77 opac- 
ity is dominant: for a cutoff energy above 1 GeV for large 
£'rad and above 300 MeV for lower Emd- For lower cutoff en- 
ergies, the Thomson opacity is now dominated by secondary 
leptons at low f^ad and primary electrons at high frad- The 
corresponding minimum Lorentz factor is not very different 
from the previous case. In panel (c) , we study the influence of 
the peak energy by using /3 = —2.5 and -Ep,o = 100 keV. The 
main effect is that the Thomson opacity is now always domi- 
nated by the contribution of primary electrons. We have also 
studied the effect of a change in the variability timescale: it 
mainly moves the limits towards higher Lorentz factors when 
the variability timescale decreases. 



3.7.4 77 opacity due to the interaction with the 
circumburst medium 

Another mechanism which could be a cause of 77 annihila- 
tion in GRBs are scatterings of 7-ray photons by electrons in 
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Figure 10. Opacity constraints on the minimum Lorentz factor in GRB outflows. The limits for the outflow to be optically 
thin for the Thomson opacity due to primary electrons (solid thick line) and secondary leptons (dashed thick line) are plotted in the 
radiated energy-Lorentz factor plane, using Eq. JSSj and Ea.l l56l l with Atvar = Is and /-, = 0.02. Note that £rad ™ the y-axis is the total 
radiated energy on a timescale Afvar- When necessary, the radiated energy above i5p,o is deduced from f^adi o and /3. Lines of constant 
cutoflf energy due to 77 annihilation are also plotted in thin line for different values of Ecut using Eq. I I54II : the dotted part of the line 
corresponds to the region where the outflow is optically thick for the Thomson opacity. Two dot-dashed lines indicate the limits where 
Ecat = itself and Ecut = Ep^Q- For Ecut < Ep^o, the lines have been extended using a low-energy slope a = —1.1. For Ecut > Egeif, we 
consider that the cutoff saturates at Eg^f = 2rmoC^. The three panels correspond to three different sets of values of /3 and Epo- 



the circumburst medium (| Thompson fc Madau 200CI'). This 
proces s has been investigated in details by iBeloborodovl 
l|2002l ). After their emission in the jet, 7-ray photons will 
progressively overtake the relativistic outflow due to a 
difi^erence ~ c/2r^ in velocity, and sweep the ambient 
medium which has not been shocked yet. Under appropriate 
conditions, this can lead to an efficient pre-acceleration of 
the circumburst medium, as photons deposit energy and 
momentum via Compton scatterings. The small fraction of 
scattered photons tends to be isotropized by the scatterings, 
which results in efficient 77 annihilation with unscattered 
photons (due to the possibility of head-on interactions), 
and therefore to the pair-enrichment of the circumburst 
medium. Both Compton scatterings and 77 annihilations 
are strongly reduced when the radius increases (dilution) 
and when the motion of the medium becomes relativistic 
(beaming). The dynamics of the deceleration (which is 
delayed) and the emission of the shocked external medium 
(which is lepton rich) are modifled by this mechanism : first 
the relativistic outflow is expanding in an empty cavity, as 
the circumburst medium is moving faster; then the external 
shock forms but is weaker due to the relativistic motion of 
the pre-accelerated material ; at larger distances there is 
no pre-acceleration any more but the circumburst medium 
remains pair-enriched; finally the standard deceleration is 
recovered. 



than a characteristic radi us R-,-,,c^t given by (see equation 
(70) of lBeloborod"wll2002l ): 



-ft'Y 



(6 ■ 10" ^ 10^") 



-Eri 



Me 



6-2/3 



(57) 

where -Brad, 53 is the radiated prompt 7-ray energy of 
the GRB in unit of 10^^ erg, ^e is the mean particle 
mass per electron in the ambient medium (/ic = 1 for 
ionized hydrogen), A is the normalization of the wind 
density profile (see Ea. (fi2)) ) and ^, = 5 • lO^^g-cm"-^ 
is the value of A expected for a typical Wolf-Rayet star wind. 

We note that the escape radius i?csc of the prompt pho- 
tons is usually larger than their emission radius Re, except 

for photons radiated at the very front of the jet. Typically, 

— 2 
-Rose — -Re + 2r Ac, where Ac is the width of the outfiow 

located in front of the emission region. Taking into account 
this effect, we compute the escape radius of all emitted pho- 
tons in the synthetic burst used to model GRB 080916C in 
Fig. [S] We obtain a minimum value of -Rcsc — 10^^ cm. For 
/3 = -2.2, -E,ad,53 = 88 (as in GRB 080916C), and pc = 1, 
we get i?T^,cxt ==; (2.2 -^ 3.6) ■ 10^^ {A/A,f-^^^ cm. There- 
fore, the condition -Rcsc ^ ^77,ext to avoid a break in the 
prompt spectrum is fulfilled as long as A/ A, < 0.2 — >■ 0.05. 



Such a mechanism may affect the prompt high-energy 
spectrum of the GRB, due to the 77 annihilations be- 
tween scatter ed and unscattered 7-ray photons. As shown by 
IBeloborodovl (j2002), when effective, these annihilations lead 
to a break at 5-50 MeV in the prompt GRB spectrum. For 
a typical high density wind-like medium, the annihilations 
rate is large enough to produce a break if the prompt GRB 
photons enter the ambient medium at a radius -Rcsc lower 



4 CONSEQUENCES OF DISTINCT EMISSION 
REGIONS FOR MEV AND GEV PHOTONS 

4.1 Are GeV and MeV photons produced in the 
same region ? 

It has been proposed in several recent studies that the 
delayed onset and/or the long-lasting tail of the high energy 
emission could be an evidence of two different regions for 
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the emission of MeV and G eV photons. An extreme version 
is t he scenario proposed bylKumar fc Barniol DuranI (|2010|) 
and iGhiseUini et alj ([201(1) where the whole GeV emission 
(prompt and long lasting) is produced by the external shock 
during the early deceleration of the relativistic outflow. 
Note that this scenario leads to strong constraints on 
the density and mag netization of the external medium 
ijPiran fc Nakan [2010) and that the observed temporal 
slope of the long-lasting high-energy emission would imply 
a strongly pair-enric hed medium as discussed in i]3.7.4l 
l|Ghisellini et al.ll2O10l ). 

Even in scenarios where the prompt GeV emission has 
an internal origin, a partially distinct emission region could 
be due to a spectral evolution of the prompt mechanism. 
For instance, in the framework of internal shocks, the evolu- 
tion of the physical conditions in the shocked region during 
the propagation of a shock wave leads to an evolving effi- 
ciency of the IC scatterings, depending on the importance 
of Klein-Nishina corrections. This naturally leads to a vari- 
able high-energy component following wit h a delay the main 
(Band) component in the MeV range (jWang et al.l l2009l : 
iBosniak et al.ll2009l : iDaigne et al.ll201lh . Successive gener- 
ations of collisions in a variable outflo w can also lead nat- 
urally to different emission regions (|Lill2010l ). An evolution 
in the microphysics of the acceleration process could also be 
responsible for some spectral evolution in scenarios where 
the re is a dominant h adronic component at high energy (see 
e.g. lAsano et al.ll2009l ). Finally, two emission regions are nat- 
urally expected in photospheric models, as it is often as- 
sumed that the main (Band) component has a photospheric 
origin and that internal shocks or magnetic dissipation oc- 
curring at larger dista nce produce an additional componen t 
at high energy (see e.g. lToma et al.ll201ll : IVurm et al.i201l[ ). 



4.2 Loosening the constraint on Fn 



As discussed in 



I Zhao et all (|2oTl); IZou et al.l (l201ll ): 
iRacusin et al.l (|201ll ). the possibility for the GeV photons 
to be produced in a different region than the MeV pho- 
tons can loosen the constraint on the minimum Lorentz 
factor in GRB outflows. To investigate this effect, we con- 
sider the same synthetic GRB as used in ^3.31 to model 
time bins 'a' and 'b' of GRB 080916C. We focus on the 
onset of the GeV component, which occurs at iobs, onset = 
iobs.trig +0.67(1 + 2) s, whcrc fobs, trig correspouds to the 
observer time of the first MeV photons. The MeV photons 
observed at fobs, onset are emitted at radius -RmcV and it is 
assumed that the emerging GeV photons observed at the 
same time were emitted by material moving with Lorentz 
factor FceV (velocity /^gevc) at radius RgcV and time fccv 
(source frame) with fccV — Rgcv/c = fobs,onset/(l + z). The 
fiash of GeV photons emitted at RgcV is assumed to have a 
power-law spectrum with photon slope /3 = —2.2. We define 
a latitude-averaged 77 opacity for GeV photons of energy 
EgbV by 



y(EGoV) 



/ 0(60)1-" sin eodOc 



(58) 

where Tj-,{EgcV,Qc) is the opacity seen photons emitted 
at colatitude Oo and D (Oc) ~ (FgoV (1 — /3gcV cosOc))~ 



is the corresponding Doppler factor. The contribution 
of each colatitude to the mean value is weighted by the 
corresponding fluence, leading to the 1 — /? exponent. 

We plot in Fig. [TT] the evolution of the latitude aver- 
aged 77 opacity r.y-y at 16 GeV (source frame) as a func- 
tion of RgcV for RgcV > RmbV and for different values 
of the Lorentz factor FgoV- When GeV and MeV photons 
are emitted at the same location, we flnd that r^-y < 1 for 
FgcV > rGcV,min, same zone ^ 340, i.e. the Same limite as in 
i]3.3l When RgcV increases, the opacity T-y^ decreases as ex- 
pected, which loosen the constraint on the minimum Lorentz 
factor rniin,GeV of the material emitting GeV photons : 



Rgcv/Rtv 



1.2 



5.1 



13 



rmin,Gev/rn 



I, GeV, same zone 



0.59 0.29 0.15 



This follows approximatively the dependency on Rc/Ro 
found in H2.2.2I i.e. r-,^ oc {R^/Rof'^'^''^'^ when R^ > Ro, 
leading to Fmin.GeV oc {Rc/Ro)~ ■ As shown in ^3.3\ the 
detailed modeling of the 77 opacity in a scenario where GeV 
and MeV photons are emitted in the same regions leads to a 
reduction of the minimum Lorentz factor by a factor ~ 2-3 
compared to single zone models. The calculation presented 
here shows in addition that the minimum Lorentz factor 
can be reduced further more by another factor ~ 2 — 8 for 
FgoV if GeV emission becomes efficient at a radius larger 
than for MeV photons. Assuming that the radiated energy 
at RgcV is not larger than the radiated energy at RmbV, we 
have checked that the outfiow remains optically thin for the 
Thomson opacity due to primary electrons and secondary 
leptons at RgcV in the case shown in Fig. 1111 

Note that the loosening of the constraint on the 
minimum Lorentz factor does not apply to models where 
GeV photons are entirely due to the external shock. Indeed, 
the small value of fobs, onset implies an early deceleration. As 
the isotropic equivalent energy of GRB 080916C is huge, 
this leads to a minimum Lorentz factor F > 10'^ in the 
outflow, which is more constraining that the 77 opacity 
limit. 

The discussion of the effect of a distinct GeV emission 
region presented here is quite simplified and some limita- 
tions should be kept in mind. If an additional GeV com- 
ponent could be firmly identified in GRB 080916C, the 
maximum energy i5MeV,max of photons associated with the 
main component should be taken into account to derive 
a new constraint Fmin on the Lorentz factor of the out- 
fiow during the MeV emission phase. We have assumed 
here F = 340 as derived in i]3.3l using -BMeV.max ~ 16 
GeV (source frame) but Fmin will be reduced by a fac- 

13 + 1 

tor (-BMeV,max/16 GeV) 20-1) if iJMeV,max is lower. There 
is one further complication: the component produced at 
RgcV extends probably in the soft gamma-ray range, as sug- 
gested by the observation of a soft excess correlated with the 
high ener gy component in some GRBs s uch as GRB 090926 
( Ackcrma mi et all l201ll'). GRB 090926 (lAbdo et all l2009ah 
and GRB 090510 (JAckermann et al.ll2010h . It has been as- 
sumed here that the annihilation rate of GeV photons with 
the seed photons produced at RgsV is negligible compared 
to the annihilation rate with MeV photons produced earlier. 
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Figure 11. Two emitting region scenario. The opacity t-/^ 
seen by 16 GeV photons (source frame) observed at tobs, onset — 
tobs,trig — 0.67 (1 + ^:) s (see text) is plotted as a function of their 
emission radius HgcV for different values of the Lorentz factor 
ToeV = 50, 100, 200, 400 and 800 from top to bottom. 



This is however not necessarily the case depending on the 
relative intensity of the two components. Clearly, a detailed 
modeling of the emitted spectrum is necessary to investigate 
such effects. This is beyond the scope of this paper and we 
leave to a forthcoming study the coupling of the formalism 
presented here to compute the 77 o pacity with a detaile d 
radiative model such as developed bv lBosniak et al.l (|2009|). 



5 CONCLUSIONS 

Photon-photon annihilation is an important process in 
GRBs. It has been known for a long time that it provides 
one of the main constraints on GRB models, i.e. the mini- 
mum Lorent z factor in the jet. In m ost GRBs, this leads to 
Tmin ~ 10^ iLithwick fc Sarill200ll ). The detection of GeV 
photons in a few GRBs by Fermj-LAT has l ed to much 
strong er limits, of the order of Fmin — 10"^ (JAbdo et al.l 
l2009bl ). However these estimates are based on highly sim- 
plified single zone formulae, which assume an isotropic ra- 
diation field in the coming frame of the outfiow. We present 
here a detailed calculation of the 77 opacity in GRB outfiows 
taking into account the exact time-dependent anisotropic 
photon field produced by multiple relativistically geomet- 
rically thin emitting regions. This calculation can be im- 
plemented numerically, for any model prescription of the 
prompt emission mechanism, as long as the emitting regions 
can be considered as expanding thin spherical shells. We 
present results obtained in the framework of the internal 
shock model, where the emitting regions are shock waves 
propagating within the outfiow. Our calculation has been 
validated by a comparison w i th the previous semi-analytical 
study made bv lGranot et al.l (|2008l ). We have also estimated 
the opacity due to Thomson scatterings by ambient electrons 



in the outflow, to pairs created by 77 annihilation and to 
77 annihilation with prompt photons back-scattered by the 
external medium. For all the examples presented in the pa- 
per, the 77 opacity is the dominant term (see i]3.7p . We have 
obtained the following results: 

(i) Spectra l shape of the 77 cutoff: as described by 
iGranot et al.l l|2008l ). due to a time evolution of the 77 
opacity, the spectral shape of the 77 cutoff in time- 
integrated spectra is not expected to be exponential, but 
rather power-law like. We illustrate this effect in a synthetic 
burst in ^3.41 This spectral shape makes more difficult the 
identification of the cutoff in fermz-LAT spectra. T he main 
candidate is GRB 090926 A (JAckermann et ahllioTI] ). where 
the number of photons detected above the cutoff at ~ 1.4 
GeV is too low to allow for a precise characterization of the 
spectral shape. Nevertheless, in most .Ffermi-LAT GRBs, 
there is no evidence for an observed cutoff, which puts 
interesting constraints on the outflow. 

(ii) Minimum Lorentz factor in GRB outfiows: when tak- 
ing into account the exact photon field, the 77 opacity for 
GeV photons is reduced compared to single zone isotropic 
estimates, leading to minimum Lorentz factors which are a 
factor of ~ 2 — 3 lower. As an illustration, we present in i]3.3l 
a synthetic burst that reproduces well the observed features 
of time bins 'a' and 'b' of GRB 080916C with a mean Lorentz 
factor of ~ 340. This minimum Lorentz factor is obtained 
assuming that there is only one component in the emitted 
spectrum, responsible both for the MeV and the GeV pho- 
tons. This is not necessarily the case as several Fermi-LAT 
burst show evidence of an additional high-energy compo- 
nent. We show in i ]4.2l that the minimum Lorentz factor can 
be reduced by another factor 2 to 8 if the GeV photons 
are emitted at larger radius than the main component. This 
study clearly illustrates the need for a detailed modeling to 
constrain the Lorentz factor in GRB outflows. When it is 
not possible, a reasonably accurate estimate of Fmin can be 
obtained from the following formula 
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(59) 



where Ci ~ 4 ■ 10~^, z and D-l{z) are the redshift and 
the luminosity distance of the source, Atvar is the ob- 
served variability timescale, RgcV / RmcV is the ratio of the 
radii where the GeV and MeV components are emitted, 
and where the high energy spectrum (over a duration 
~ Aivar) is assumed to follow a power-law with photon 
index /3 above an observed characteristic energy Ec : 
F{E) = F{Ec){E/E,f (ph.cm-^keV-i). Energy Smax 
is the observed energy of the most energetic detected 
photons. The integral I{/3) is defined by Eq.(1221l in ^^^ 
As usually the spectrum is measured over a time interval 
Atspoc which is larger than the variability timescale Atvar, 
the normalization F{Ec) entering in Eq. l|59p (fiuence at 
energy Ec in ph.cm~^.keV~^) must be corrected by a factor 
F{Ec) = F{Ec) X (Aivar/Atspcc). This equation can be 
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directly applied to ■Fermi-L AT observations and generalizes 
the usual formula given by lAbdo et al.1 l|2009bl ) by intro- 
ducing two corrections: (1) a more accurate normalization 
taking into account the anisotropy of the radiation field and 
including a numerical factor Ci obtained from the compar- 
ison with numerical simulations in i i3.3l (2) the possibility 
to take into account two different emitting regions for MeV 
and GeV photons. The standard limit is obtained with 
RgcV / RmbV = 1 (same region) : then the denominator in 
Ea. (|59p equals 1. The radius RmcV is estimated from the 
variability timescale by 7?MeV — r^cAivar/(l + z), which is 
valid for most models of the prompt emission. The radius 
RceV is difficult to constrain without a detailed model of 
the high-energy emission mechanism. If GeV photons have 
an internal origin, an upper limit for RgcV is given by the 
deceleration radius, fn the future, a measurement of the 
variability timescale in the GeV lightcurve could provide a 
better estimate of this radius. 

(iii) Delayed onset of the GeV emission: due to the 
variable nature of the GRB phenomenon, the 77 opacity 
is expected to be strongly time dependent. This can lead 
naturally to a delayed onset of the high-energy component, 
if it is initially highly absorbed. For instance, in internal 
shocks, the emission radius increases during the propagation 
of a given shock wave, which favors such an evolution. We 
show several examples of synthetic bursts where a delayed 
onset of the GeV component is observed, fn the case of 
internal shocks, the delay before the onset is comparable 
with the observed duration associated with the propagation 
of the shock wave that produces the seed photons for 77 
annihilation. This duration has to be distinguished from 
the shortest timescale of variability, fn the synthetic burst 
which models time bins 'a' and 'b' of GRB 080916C, we 
obtain a GeV onset delayed by ~ 5 s, comparable with the 
observed value, while reproducing variability on shorter 
timescales (~ 0.5 s) in the MeV lightcurve. A time evolving 
77 opacity appears as a good candidate to explain the 
observed delayed onset of the GeV emission. Note that this 
effect is obtained here assuming a single component in the 
emitted spectrum. More complex spectral e volution, such 
as a varying inverse Compton component IjBosniak et al.l 
I2OO9I ) , may be an additional source of delay. 

(iv) Smoothing the variability in the GeV lightcurve: in 
highly variable outflows, the short timescale variability in 
the GRB lightcurve is expected to be associated with emis- 
sion regions at low radius. Therefore, the associated high- 
energy emission may be highly absorbed by 77 annihilation 
due to a denser photon field. This predicts a smoothing of 
the lightcurve in the GeV range. We illustrate this effect in 
the framework of the internal shock model in H3.6I For this 
reason, the variability of the GeV lightcurve is an indica- 
tor to use carefully if one wants to distinguish between an 
internal and an external origin for the GeV emission. 

Our results loosen the constraints on models of GRB cen- 
tral engines as Lorentz factor of 10'' do not seem to be 
required, even in the most extreme bursts observed by 
Fermi-ljAT. In addition, they solve another contradiction. 
For Lorentz factors of 10"^, the deceleration by the exter- 
nal medium occurs very early. Then, it becomes difficult 



to interpret Fermi-LAT GRBs such as GRB 080916C ff 
the pro mpt emission is radiat ed well above the photosphere 
(see also lZhang fc Pe'erll2009l '). On the other hand, the early 
steep decay observed in X-ray afterglows by Swift-XRT is 
a strong eviden ce that the prompt emission phase ends at 
a larpe radius ([Lyutikovl l2006l: iLazzati fc Begelmaa l2006l : 



iKumar et ai]|2007l : lGenet fc Granotll2009l ). With lower lim- 
its on the Lorentz factor as obtained in the present paper, 
bright Fermi-LAT GRBs such as GRB 080916C are con- 
sistent with models, including internal shocks, where the 
prompt emission is produced between the photospheric and 
the deceleration radii. 



APPENDIX A: PHOTON FIELD CREATED BY 
A SPHERICAL FLASH 

Let us consider the interaction between a spherical flash and 
a high energy photon (the geometry is illustrated in Fig. [U 
see also the caption of the figure for the adopted notations). 
The emissivity in the source frame is given : 



where the Doppler factor T> is defined by 
V - 



V' r (1-/3 cos 5) 



(Al) 



(A2) 



The specific intensity 7,y [erg-cm~^-s~^ •Hz"'^ -sr"^] created 
by the spherical flash at the interaction point I (radius Ri 
and time ti) in direction a is given by 



Iv {Ri,ti,a) = / > s,ti 



ds , 



(A3) 



assuming negligible absorption at frequency v. Indeed, high 
energy photons interact preferentially with seed low energy 
photons close to the threshold with -0 ~ l/F, i.e. 



^ ^"^ ^ 1.04 (^ 
10 MeV 1 GeV VlOO 



(A4) 



As the number of emitted photons at energy E, N{E), de- 
creases as a power-law with energy (typically N{E) oc E^ 
with /3 ~ -2 to -3 above 1 MeV), the ratio N{EiiE)/N{E) 
is always small and the fraction of absorbed seed photons 
negligible. The radius r is related to s and a by 



r — Ri + S ~2Riscosa, 



which leads to 



(5 (r — _Ro) ds ■ 



Ro 



s — Ri cos a 
The same relation applied at r = _Ro (S ~ s) leads to 
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Therefore, from Eq.® and Ea. (|A3|l . we obtain 
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As expected from the flash geometry, the intensity is positive 
for a unique direction a = qi. Then, the photon density 
[ph ■ cm~ • keV~ ■ sr~^] is given by 



nn{E) = 



he hi/ 



(A9) 



the following approximate formula: 
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APPENDIX B: APPROXIMATE FORMULA 
FOR THE 77 OPACITY 

In this appendix, we derive the approximate formula given 
by Eg. H'iSp for the 77 opacity seen by a on-axis high en- 
ergy photon emitted at 7?e interacting with a spherical flash 
radiated at _Ro and having a power-law spectrum. Using 
-Re = i?o + /3oc (tc —to) and introducing normalized coor- 
dinates X = (_Re — Rti) /Ro et X = £/Ro, the opacity given 
by Eq.((23l) reads 

r,,(£HE) = X(/3)ro^o f y^^l rj+^ 



EueEp.o 
1 X 



dx 



X 1+X + x 



+ x 



V^-^ (1-cosV)"'' 



(Bl) 



[(l-fX)(l+f 



f)]" 



(B2) 

which leads to Ea. H25p . This approximate formula is com- 
pared to the exact calculation in the top- left panel of Fig. [S] 
The behaviour of the opacity with the emission radius of 
the high-energy photon can be simply understood from the 
Taylor expansions above : the integral is dominated by the 
product of the dilution, Doppler et interaction angle terms. 
When the high energy photon is emitted just after the flash 
(-Ro — >■ -Ro, i-e. X — > 0), the dilution term is negligible, the 
Doppler term is optimal as 5 ~ I/Fq and the interaction 
angle term is also optimal as i/; ~ I/Fq. Therefore, 77 an- 
nihilation is efficient. At larger emission radii (_Rc 3> -Ro, 
i.e. X — >■ -l-oo), the dilution factor (oc X^^) reduces the 
seed radiation field, the Doppler term (oc X^~^) also re- 
duces strongly the seed radiation field as 5 becomes large 
compared to 1/Fo, and the interaction angle term (oc X^) 
strongly reduces the probablity of interaction as t/j — > 0. This 
explains the steep decay of Tjj oc X^^^~^' for X > 1 . 



Several contributions appear in the integral 



• Dilution term : the first factor, {1 + X + x)~ , ex- 
presses the dilution of the radiation from the flash since its 
emission ; 

• Propagation term : the second factor, {X/Po + x)^ , 
takes into account the propagation of the high energy photon 
up to location of the interaction ; 

• Doppler term : the third factor, T>^~^ ^ is due to the 
Doppler effect applied to the seed photons. 

• Interaction angle term : the fourth term, (1 — cosi/;)" , 
expresses the effect of the interaction angle ip on the effective 
cross section. 

The product of these four terms is a strictly decreasing 
function, so that the integral is dominated by the con- 
tributions at low X. For x <^ X and Fo 2> 1, we have 
{l + X + x)~^ ~ (l + X)-' (dilution), X/iX/Po+x) ~ 
/3o ~ 1 (propagation), D^-" ~ (2Fo/(2 -F X))^"" (Doppler) 
and (1— cos?/))~ ~ (2Fo(l + X)) (interaction angle), 
where we have used the Taylor expansion of the Doppler 
angle 5 and the interaction angle t/j : 



5 = 

lb ■^ 


VI + x 

Fo ' 
1 




Vi + ^ro 



The product of the four terms is therefore approximatively 
constant for x <^ X, before starting a steep decay. Keeping 
only the contribution for ^ a; < X to the integral gives 
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